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ABSTRACT 

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 
125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical 
continuum sampling ever obtained, with a mean sampling rate <0.5-day. Approximately daily HST 
UV sampling was also obtained. The UV/optical light curves show strong correlations (umax = 
0.57 — 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a 

T OC A4/3 

wavelength dependence, with a normalization that indicates an unexpectedly large disk 
radius of ^ 0.35 ± 0.05 It-day at 1367 A, assuming a simple face-on model. The [/-band shows a 
marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer 
continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray 
correlation is weaker (rmax < 0.45) and less consistent over time. This indicates that while Swift is 
beginning to measure UV/optical lags in general agreement with accretion disk theory (although the 
derived size is larger than predicted), the relationship with X-ray variability is less well understood. 
Combining this accretion disk size estimate with those from quasar microlensing studies suggests that 
AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses. 
Subject headings: galaxies: active - galaxies: individual (NGC 5548) - galaxies: nuclei - galaxies: 

Seyfert 
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1. INTRODUCTION 

Because of their great distances and small sizes, the 
central regions of active galactic nuclei (AGN) can¬ 
not be resolved directly with current technology. Thus 
it is necessary to use indirect methods to gain infor¬ 
mation about AGN structure and physical conditions. 
Variability studies, along with gravitational mic r olens- 
ing (e.g.. [Morgan et a .11201(1 iM osauera et al l [ 201 .^ 


iBlackburne et al.l 12014 . I.Timinez- Vicente et al.l [201411 . 

have emerged as powerful techniques for probing the cen¬ 
tral regions of AGN. 

In p articular, the “reverberatio n mapping” (RM) tech¬ 
nique (|Blandford fc McKeelll982[ l has proven quite effec¬ 
tive at taking advantage of strong AGN line and contin¬ 
uum variability to probe the structure of the broad emis¬ 
sion line region (BLR). The fundamental idea of RM is 
that if the variability in band B is powered by variabil¬ 
ity in band A, with only light travel times affecting the 
light curves, then variations in band A will be seen in 
band B, but delayed and smoothed by the size and ge¬ 
ometry of the latter emitting region. The first unambigu¬ 
ous application of RM came in an lUE campaign on the 
Seyfert 1 galaxy NGG 5548, which found that variations 
in the driving UV continuum (band A in this picture) 
were highly correlated with those in emission lines such 
as C IV line (band B). The line variations lagged the con¬ 
tinuum by ^10 days, indicating tha t the G iv-emittin g 
region was of order 10 It-days in size (jGlavel et al.lllOOll l. 
Optical emission lines showed similarly strong correlation 
but w ith larger lags. For ex ample H/I showed a lag of ^20 
days (|Peterson et al.l 1199111 . indicating a stratified BLR 
in which higher-ionization lines are formed closer to the 
central engine. The distance estimate, when combined 
with the line width, allows estimation of the mass of the 
central black hole. For NGC 5548 t he current best mass 
estimate is Mrh ^3 .2 x 10^ Mq (iDennev et al. I |2M^ 
iPancoast et al.ir2014f) . This technique is now a standard 
tool for AGN astronomy, yielding BLR size, stratification 
information and black hole mass est imates and physical 
conditions for ^50 AGN (see, e.g., iBentz fc Katzl[2015l 
for a recent compilation). For a more extensive general 
discus sion of BLR RM, plea se see the first paper in this 
series (|De Rosa et al.l[2015t Paper I hereafter). 

The structure and physics of the central engine that 
produces the continuum emission is currently less well 
understood than the reverberation-mapped BLR. Ob¬ 
servations and accretion disk theory both suggest that 
the inner accretion disk/corona region emits short wave¬ 
length continuum (X-ray, ultraviolet [UV hereafter] and 
much of the optical), which then illuminates and ionizes 
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the gas in the more distant BLR and beyond. The pre¬ 
vailing picture is that the black hole is surrounded by a 
small, hot (T ~ 10® K) and relatively spherical corona 
and a larger, cooler (T /nay ^ 5 x 10^ K) and re latively 
flat accretion disk fe.g.. lHaardt &: MaraschillQQlll . Grav¬ 
itational leasing studies also indicate that this putative 
corona is small enough t o be considered point-like rela¬ 
tive t o the disk (^ 5Rc|„ uilDai et al. 11201 3iJMorgane^jl 


_ nhiJiJaijerajjizuiijjjviorgane^ 

2 OI 2 L iMosouera et al.T 1^1,'U IBlackburne et al.l I20T_ 
Blackburne et al.l 1201511 . The energy released by the 


accretion process heats both the optically thick disk 
- producing the thermal UV/optical emission - and 
the corona, which in turn can illuminate the disk as 
an external heating source. The fraction of the en¬ 
ergy that goes into heating the corona has not been 
established and therefore it is not clear whether the 
disk is mainly heated internally or externally. In ei¬ 
ther case, however, the disk is expected to have a strati¬ 
fied temperature structure with the hotter, UV-emitting 
regions closer in and the cooler, optically-emitting re¬ 
gions farther out. Quasar microlensing studies find 
that accretion disk size s increase with wavelength (e.g., 
iPoindexter et al.l 1200811 . supporting this picture. We 
note however that this picture contains important un¬ 
reconciled discrepancies. For instance gravitational leas¬ 
ing disk sizes are typically r eported to be a fact or of 
^4 larger than predicted le.g.. iMorgan et aDl2010l l. and 
the observed UV spectrum is too steep with a Lyman 
discontinuity that is typically smaller than predict ed or 
not seen (e.g. iKoratkar &: Biaejll999L lGollinll2001ll . Re- 
cent improvements in AGN accretion disk models (e.g. 


IDexter Kr. Ae'olll2011 

1 may 

see also 1 Antonuccil ( 

20131) 


Just as RM of the BLR allows us to estimate the 
distances at which each line is formed, RM of the ac¬ 
cretion disk could allow us to constrain the temper¬ 
ature struct ure of the disk and test the standard a- 
disk model (jShakura fc Sunva^ll973l or any other pre¬ 
dictive model). Repeated efforts have been made to 
implement RM of the accretion disk by correlating X- 
ray light curves gathered with space-based observato- 


observatories ( 

e.ff.. lEdelson et al 

119961 INandra et all 

19981 ISuganu 

ma et al. 1 120061 

Arevalo et al. 2008. 

Arevalo et al. 1 

20091 IBreedt et al.l 120091 IBreedt et al.l 

20101 IGameron et al.l 2012L IGliozzi et al.||2013|). How- 


ever the optical time resolution of these early experi¬ 
ments was typically limited by the diurnal cycle to AT > 
1 day, resulting in lag measurements that were sugges¬ 
tive but not statistically significant (1 — 2a), although 
often in the expected sense, with X-rays leading the op- 
tical. Other experim ents from this peri od used HST 
(lEdelson et al.l 1200011 or XMM-Newton (|Mason et al.l 
1200211 to attain finer optical time resolution at the 
cost of shorter monitoring periods (~l-2 days). Again 
the results were suggestive but inadequate to make a 
definitive lag measurement. Further, these experiments 
typically only sampled a single optical or UV band, 
and thus were unable to explore temperature stratifi¬ 
cation in the disk. Ground-based multicolor optical/IR 
studies also yielded tentative evidence of shorter wave¬ 
length variatio ns leading longer wavelength variations 
in some AGN (j Sergeev et al.l 20051 iGackett et al.ll2?)f)7l 




















































































AGN STORM. 11. Swift Observations 


3 


iLira et aDl2011h . 

The unique capa bilities of the Swift observatory 
(|Gehrels et al.l[200ll . originally optimized to detect 7 - 
ray burst counterparts, are also ideally suited for AGN 
monitoring. Its rapid slew/acquisition times and large 
sky coverage make it feasible to sample AGN light curves 
(which show variability over a broad range of tempo¬ 
ral frequencies) with high cadence over a long dura¬ 
tion. Further , the coaligned Ultr aViolet-Optical Tele¬ 
scope lUVOT:iR oming et al. 11200^ and X-Ray Telescope 
(XRT : lBurrows et al.ll200!Tfl cover the entire energy range 
of interest (the X-ray/UV/optical) with a single space- 
based telescope, so data quality is no longer limited by 
the diurnal cycle or weather. 

This is leading to important advances in accretion disk 
RM, as highlighted by the success of two recent Swift 
AGN monitoring campaigns. After detecting from the 
ground that the relatively no rmal galaxy NGC 261 7 had 
transitioned into a Seyfert l. lShaopee et all (1201411 used 
Swift to cover a ^50 day period with approximately daily 
cadence, generally in all six UVOT filters. Ground-based 
optical/infrared coverage of the first part of this period 
was obtained at a lower cadence. iMcHardv et al. I (|2014f l 
analyzed 359 “visits” (separated by 1 orbit or longer; 
see Section 2.1) to the archetypical Seyfert 1 galaxy 
NGG 5548, over -2 years (2012 Feb - 2014 Feb). Ap¬ 
proximately 20% of the visits utilized all six filters. In 
both cases, the data show significant (> 3cr) interband 
lags throughout the UV/optical, with increasing lags to 
longer wavelengths, consistent with a dependence 
as predicted by the standard a-disk model under the as¬ 
sumption that time lags are dominated by light travel 
times. 

The experiment detailed herein combines Swift’s pow¬ 
erful capabilities with simultaneous, intensive UV spec¬ 
troscopic monitoring by HST to yield the densest X- 
ray/UV/optical coverage - in both time and wavelength 
- ever obtained. NGC 5548 is the target of this cam¬ 
paign. Figure 1 shows that this campaign yields a factor 
of ^2.5-4 improvement in the number of UVOT filter 
data pairs available for correlation compared to the two 
best previous campaigns. This provides superior power 
to measure small (< 2 -day) interband lags with high pre¬ 
cision. 

The result, presented in this paper, is a clear mea¬ 
surement of lags across the entire UV/optical range, 
with shorter wavelength bands leading the longer wave¬ 
length bands. The timescales generally increase to longer 
wavelengths as expected for a standard a-disk, but di¬ 
rect fitting indicates a larger than expected disk (^ 
0.35 ± 0.05 It-day at 1367 A). The [/-band lag is slightly 
longer than expected from the fits, apparently consis¬ 
tent with conta mination from BLR cont inuum emission 
as predicted by iKorista fc GoadI (|2001ll . The X-rays 
show a relatively weak and less coherent relation to the 
UV/optical. Finally combining RM and microlensing 
disk size estimates suggest that disk size scales roughly 
linearly with black hole mass over a wide range of masses. 

This paper is organized as follows. Section 2 describes 
the observations and data reduction. Section 3 presents 
cross-correlation analyses applied to these data. Section 
4 discusses the theoretical implications of these results, 
and Section 5 concludes with a brief summary of this 


work and implications for the future. 

2. OBSERVATIONS AND DATA REDUCTION 
2 . 1 . Observations 

The t arget of this experiment, NGC 5548 {z = 
0.01717: Ide Vaucouleurs et al.lll99l]l . shows strong, reli¬ 
able variability across the entire X-ray/UV/optical wave¬ 
length range accessible to Swift. It is also among the 
brightest AGN in the sky at these wavelengths. In 2014 
February - June, Swift executed a monitoring campaign 
on NGC 5548 that was ground-breaking in two respects: 
1) it was comprised of 360 separate visits over a ~4 
month period, of which 282 successful visits were ob¬ 
tained, for a sampling rate (after removing bad data) bet¬ 
ter than one vi sit every 0.5-day , and 2) it utilized all six 
UVOT filters (iPoole et al.ll2008ll in each visit, with 239 
(84%) providing usable measurements in all six filters. 
(For the purposes of this paper, a visit is defined as an 
observation in which at least one UVOT filter measure¬ 
ment is obtained. Multiple observations within a single 
^96 min orbit are combined to form a single visit.) This 
entailed a significant commitment of spacecraft resources 
given the limit of 500 time-critical non-GRB guest inves¬ 
tigator visits per yea r and the desire to minimize wear 
on the filter wheelo 

In addition, a parallel HST emission-line RM campaign 
yielded daily UV spectroscopic monitoring of NGC 5548 
over a slightly longer period (see Paper I). This provided 
mutual synergies: the HST 1367 A continuum light curve 
was used in the cross-correlation functions (CCFs) re¬ 
ported herein, while the Swift optical, ultraviolet and 
X-ray light curves can be used to better define the con¬ 
tinuum variability characteristics needed to understand 
the emission-line RM results. 

These observations are summarized in Table 1. Start 
and stop times for Swift observations are originally 
recorded in MET (Mission Elapsed Time; seconds since 
the start of 2001 ) and corrected for the drift of the on¬ 
board Swift clock and leap-seconds. These times were av¬ 
eraged and converted to Heliocentric Julian Date (HJD), 
the standard for this observing campaign. Through¬ 
out this paper we utilize the truncated HJD, defined as 
THJD = HJD - 2,456,000. We reduced all Swift data 
on NGC 5548 for both the UVOT and XRT, but re¬ 
stricted scientific analysis to observations taken dur¬ 
ing the intensive monitoring period, from THJD 706 
to THJD 831 (approximately 2014 Feb 17.5 - Jun 22.5 
UTC). 

The HST 1367 A data reduction is detailed in Paper I. 
The following two subsections will describe the reduction 
of the Swift UVOT and XRT data. 

2 . 2 . UVOT data reduction 

Swift observed NGC 5548 for a total of 2935 expo¬ 
sures in six UVOT filters from the beginning of the mis¬ 
sion through THJD 876. All UVOT data were repro- 
cess ed for uniformity , applying standard FTDDLS utili¬ 
ties (lBlackburnlll995t from version 6.15.1 of HEASOFtF?). 
The astrometry of each field was refined using up to 

http://swift.gsfc.nasa.gov/proposals/tech_appd/swiftta_vll 
/node42.html 

http://heasarc.gsfc.nasa.gov/ftools/ 
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TABLE 1 

Monitoring information 


TABLE 2 
UVOT DATA 


(1) 

(2) 

(3) 

(4) 

(5) 

(1) 

(2) 

(3) 

(4) 


(Jentrai 

Wavelength 

Number 

Sampling 

HJD 

Filter 

Flux 

Error 

Band 

A (A) 

range (A) 

of points 

Rate (day) 










2454270.833 

UVW2 

0.625 

0.018 

HX 

4.4 

1.2 - 15.5 

272 

0.46 

2454270.905 

UVW2 

0.632 

0.015 

SX 

25.3 

15.5 - 41.3 

272 

0.46 

2454276.539 

UVW2 

0.671 

0.017 

HST 

1367 

1364.5 - 1369.5 

121 

1.03 

2454276.606 

UVW2 

0.663 

0.017 

UVW2 

1928 

1650 - 2250 

262 

0.47 

2454283.359 

UVW2 

0.701 

0.019 

UVM2 

2246 

2000 - 2500 

254 

0.49 

2454283.427 

UVW2 

0.701 

0.019 

UVWl 

2600 

2250 - 2950 

266 

0.47 

2454283.492 

UVW2 

0.707 

0.020 

U 

3465 

3050 - 3900 

266 

0.47 

2454290.329 

UVW2 

0.802 

0.023 

B 

4392 

3900 - 4900 

265 

0.47 

2454290.378 

UVW2 

0.796 

0.020 

V 

5468 

5050 - 5800 

258 

0.48 

2454290.445 

UVW2 

0.796 

0.020 


Note. — Column 1: Observing band name. Column 2: 
Central wavelength of that band. Colu mn 3: FWHM wave - 
length range of that band, estimated from IPoole et aEl 1I2008IL 
Column 4: Total number of good data points in that band. 
Column 5: Mean sampling rate in that band. 


Note. — Column 1: Heliocentric Ju¬ 
lian Date. Column 2: Observing Filter. 
Column 3: Measured flux in units of 
10“^^ erg cm“^ A“^. Column 4: 
Measured Icr error in the same units. 
Note that this table includes all usable 


35 isolated field stars drawn from the HST GSC 2.3.2 
(|Lasker et al.ll2008h and Tycho-2 (|He)g et al.ll2000D cat¬ 
alogs, yielding residual offsets that were typically ~0.3 
arcsec. Fluxes were measured using a 5 arcsec circular 
aperture and concentric 40-90 arcsec regions were used 
to measure the sky background level. The final values in¬ 
clude corrections for aperture losses, coincidence losses, 
and variation in the detector sensitivity across the image 
plane. The galaxy contributes a fraction of the observed 
flux within the UVOT apertures (see Section 4.4) but 
no attempt was made to remove the contribution of host 
galaxy flux, as this contamination is constant and will 
not affect measurement of interband temporal correla¬ 
tions or absolute variability amplitudes. 

We screened the data to eliminate exposures affected 
by significant tracking errors. To identify observations 
with distortions in the wings of their point spread func¬ 
tion (PSF), we measured the ratio of counts in annuli 
from 5-7 and 7-10 arcsec, determined the distribution of 
these ratios for each filter and discarded any observations 
that were found to be outliers by at least 3.5o- (defined 
iteratively). In addition, we measured the PSFs of the 
isolated field stars used for astrometric refinement, flag¬ 
ging any observations for which either the average PSF 
Full width at half maximum (FWHM) differed by more 
than 1.0 arcsec from the nominal UVOT FWHM (2.2-- 2.9 
arcsec depending upon the filter. IBreeveld et al.ll2010ll or 
the average FWHM of the stellar PSF projections along 
the X and Y axes differed by more than 0.75 arcsec. All 
flagged observations were manually inspected, leading to 
the rejection of one additional exposure in which the stars 
were streaks 15 arcsec long. In total, 30 exposures are 
rejected. 

The resulting light curves exhibited occasional, anoma¬ 
lously low points, especially in the UV. Subsequent in¬ 
vestigation found that these “dropouts” occur when the 
source falls within specific regions of the detector. Data 
potentially affected by these suspect detector regions are 
identified and removed using a new methodology dis¬ 
cussed in the Appendix, eliminating 7.4% of the expo¬ 
sures. Finally we combined fluxes and errors in quadra¬ 
ture so there is no more than one data point per filter 
per orbit for any orbit in which multiple measurements 


Swift observations of NGC 5548, not 
just those from the intensive monitoring 
period. The data are sorted first by Al¬ 
ter, then by HJD. Only a portion of this 
table is shown here to demonstrate its 
form and content. A machine-readable 
version of the full table is available on¬ 
line. 

were made in the same filter. The final light curves are 
presented in Figure 2 and data from the full mission are 
given in Table 2. 

2.3. XRT data reduction 

The Swift XRT data were gathered in photon count¬ 
ing (PC]_^io^ and a nalyzed using the tools described 
by I Evans et al.l (|2009li^ to produce light curves which 
are fully corrected for instrumental effects such as pile 
up, dead regions on the CCD and vignetting. We gener¬ 
ated soft X-ray (SX; 0.3-0.8 keV) and hard X-ray (HX; 
0.8-10 keV) light curves. We utilized “snapshot” bin¬ 
ning, which produces one bin for each spacecraft orbit. 
As with the UVOT data, we averaged multiple ObsIDs 
within a single orbit in quadrature. We investigated the 
use of other bands by subdividing HX into 0.8-2.8 keV 
and 2.8-10 keV but made no change after finding the cor¬ 
relation properties of the sub-bands to be very similar to 
the original choice. 

The gap in the X-ray light curves during THJD 812-819 
(2014 Jun 4-10; Figure 2) corresponds t o the time that 
the Swift XRT was in an anomaly state (IBurrow s et al.l 

l2014aLlBurrows et al.ll2014RlKennea k. Burrowsl2014bll . 

during which time XRT was either disabled or collected 
data in a non-standard, not-fully calibrated mode. We 
excluded the data taken during this time interval from 
our analysis. We additionally excluded all visits where 
the total good integration time was less than 120 sec. 
This resulted in a final light curve having 272 X-ray 
points over the 125-day intensive monitoring period (see 
Table 1). The complete NGC 5548 XRT data are pre¬ 
sented in Table 3. 

2.4. Light curves 

http://www.swift.ac.uk/user_objects. 
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TABLE 3 
XRT DATA 


(1) 

HJD 

(2) 

HX Flux 

(3) 

HX Error 

(4) 

SX Flux 

(5) 

SX Error 

2453468.872 

0.324 

0.039 

0.093 

0.023 

2453469.005 

0.361 

0.051 

0.061 

0.021 

2453469.139 

0.282 

0.049 

0.156 

0.042 

2453470.283 

0.425 

0.057 

0.114 

0.030 

2453470.349 

0.343 

0.037 

0.075 

0.019 

2453470.418 

0.360 

0.028 

0.076 

0.013 

2453470.814 

0.353 

0.024 

0.088 

0.012 

2453470.882 

0.469 

0.043 

0.088 

0.019 

2453473.227 

0.496 

0.055 

0.175 

0.033 

2453475.169 

0.384 

0.049 

0.073 

0.021 


Note. — Column 1: Heliocentric Julian Date. Columns 2 
and 3: Measured HX flux and Icr error, in ct sec“^. Columns 
4 and 5: Measured SX flux and la error, in ct sec“^. Note 
that this table includes all usable Swift observations of NGC 
5548, not just those from the intensive monitoring period, 
sorted by HJD. Only a portion of this table is shown here 
to demonstrate its form and content. A machine-readable 
version of the full table is available online. 

Although the Swift data for NGC 5548 span many 
years, Figure 2 and Table 1 cover only the ~ 125-day 
intensive monitoring period THJD 706-831. The light 
curves are presented in order of descending frequency 
with the highest frequency band at the top and the lowest 
at the bottom. The HST light curve plays a critical role 
as the only data set not gathered by Swift. This means 
that CCFs relative t o this band will not suffe r from “cor¬ 
related errors” fsee lEdelson fc Kroliklll98^ . The HST 
light curve also has much higher signal-to-noise ratios 
and better exclusion of BLR emission than the Swift 
data, but with less than half the sampling cadence. 

3. INTERBAND CORRELATION AND VARIABILITY 
ANALYSES 

In this section we estimate the interband correlation 
and lag between continuum bands. Before performing 
these correlation analyses we detrended the data by sub¬ 
tracting a 30-day boxcar running mean. This was done to 
remove long-term trends that could potentially degrade 
our ability to measure the expected small lags. 

In all correlation analyses we reference the correlation 
of one band (the HST band) relative to all other bands 
(the eight Swift UVOT and XRT bands), restricting our 
analysis to just the data shown in Figure 2. We used the 
interpolate d cross-correlation fun ction (ICCF) as imple¬ 
mented bv iPeterson et al.1 (|2004ll . to measure and char¬ 
acterize temporal correlations within these data. These 
results are shown in Table 4 and discussed in the fol¬ 
lowing subsection. The second subsection reports the 
result of fits to these data, and the third describes our 
characterization of the variable UVOT spectral energy 
distribution (SED). 

3.1. Correlation analysis 

The traditional CCF le.g.. I.Tenkins fc Wattslll96^ re¬ 
quires evenly sampled data, but most astronomical data 
are not evenly sampled. The ICCF performs a piecewise 
linear interpolation in the reference (HST) band with a 
user-defined interpolation step of 0.1 day, and then mea¬ 
sures the correlation relative to the non-interpolated data 


TABLE 4 

Interband correlation 
COEFFICIENTS AND LAGS 


(1) 

Band 

(2) 

‘Tmax 

(3) 

Lag (days) 

HX 

0.35 

-0.66 ± 0.46 

SX 

0.44 

-1-0.08 ± 0.52 

HST 

1.00 

-1-0.00 ± 0.25 

UVW2 

0.90 

-1-0.40 ± 0.17 

UVM2 

0.87 

-1-0.35 ± 0.16 

UVWl 

0.85 

-1-0.61 ± 0.20 

U 

0.81 

-1-1.35 ± 0.24 

B 

0.74 

+1.23 ± 0.29 

V 

0.57 

+1.56 ± 0.50 


Note. — Column 1: Band 
for which correlation was mea¬ 
sured relative to HST 1367 A. 

Column 2: Maximum corre¬ 
lation coefficient. Column 3: 

Measured centroid lag and asso¬ 
ciated la error in days. 

in the other band. These data are then shifted and cor¬ 
related to build up the correlation function. In this case, 
the HST data, with an initial cadence of ^1.1 day, are 
resampled to a grid with 0.1 day spacing, and then the 
CCF of each Swift light curve is measured relative to the 
LIST light curve. 

The results are shown in Figure 3a. The third panel 
shows the auto-correlation function (ACF) of the HST 
data; all others are CCFs measured relative to the HST 
light curve, so a positive lag indicates that variations in 
that band lag behind the HST light curve. Two points 
are clear. First there is a tendency for peak lags to in¬ 
crease with wavelength: the hard X-ray band shows a 
negative lag relative to HST, the soft X-rays show ap¬ 
proximately zero lag, the lags are positive and small 
within the Swift UV (longer wavelengths lag HST) and 
the lags are positive and larger between HST and the 
Swift optical. Second the strength of the correlation is 
larger between the HST and UVOT data (peak correla¬ 
tion coefficients rmax = 0.57 — 0.90) than between the 
HST and XRT data (r^ax < 0.45). 

In order to quantify the uncertainties on the interband 
lag estimates, we utilized the flux random ization/random 
subset selection (FR/ RSS) technique of IPeterson et al. I 
(1199811 as modified bv IPeterson et al.l (1200411 to produce 
the cross-correlation centroid distribution (CCCD), as 
shown in Figure 3b. FR/RSS is a mo del-independent 
Monte Carlo technique that attempts to deal with both 
flux uncertainties in individual measurements and un¬ 
certainties due to sampling of the time series. In “ran¬ 
dom subset selection” (RSS), for a light curve of N data 
points, one randomly selects N data point without regard 
to whether a data point has been previously selected or 
not. Thus, approximately 1/e of the original points in 
the light curve are not selected in a given realization, 
and the remaining points are selected one or more times. 
For data points selected n times in a given realization, 
the uncertainty associated with the data point is reduced 
by . “Flux randomization” (FR) consists of alter¬ 

ing the observed flux by random Gaussian deviates whose 
standard deviation is equal to the flux uncertainty on the 
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data point. The CCCD is built by combining the results 
from 2000 realizations, with results that are summarized 
in Table 4. 

A number of factors can contribute to the widths of 
the histograms and thus the error estimates on the inter¬ 
band lags. These include the sampling cadence and finite 
duration of the campaign, measurement errors, and de¬ 
viations from the stationarity assumption implicit in the 
FR/RSS method. An example of the second contribu¬ 
tion could be the appearance of somewhat different lags 
at different epochs, as may be occurring with the BLR 
(see Paper I). At present it is not possible to determine 
the relative contribution of each effect. 

3.2. Lag-wavelength fits 

As discussed in the introduction, the standard model 
predicts a relationship between lag and wavelength be¬ 
cause the disk is expected to be hotter at smaller, inner 
radii a nd cooler at larg e r, out er ra dii. Following the ana l- 
ysis of lMcHardv et al.l (I2014D and lShaopee et alTl (I2014D . 
Figure 4 presents the CCF lag (r) results as a function 
of wavelength (A). We fit the wavelength dependence 
of the lags with the function t = AB{{X/Xo)^ — 1 ). 
The top three sets of panels show the effect of restrict¬ 
ing the fitting function by first setting (7 = 4/3 and 
then setting A = 0. This yields only a slight increase 
in xIj which is acceptable in all cases. For instance the 
third panel has reduced of xl = 0.98, correspond¬ 
ing to a probability value p = 0.45. Thus we conclude 
these data are fully consistent with a single-parameter 
fit, r = B{{X/Xoy^^ — 1). The fit parameter B gives an 
estimate of the size of the disk at Aq = 1367 A, the HST 
reference wavelength, assuming a face-on geometry. This 
is the most important result of this paper, discussed in 
detail in Section 4.1. 

The bottom three sets of panels explore the effect of 
excluding particular bands from the fit. The third panel 
shows the effect of excluding the HST ACF lag, which of 
course should be identically zero. This has no effect on 
the number of degrees of freedom as the fit parameter A 
is dropped as well (as discussed above). 

The fourth panel shows the additional effect of drop¬ 
ping the two X-ray lags, HX and SX. This has essentially 
no effect on fit quality, which is not surprising because 
the correlation coefficients are low and the lags have by 
far the largest errors of any waveband. That is, the X- 
ray variations do not show a clear, consistent relation 
to the UV/optical variations. The implications of this 
divergence is explored in detail in Section 4.4. 

Another new result of this experiment is that the U- 
band lag is consistently larger than predicted by the fits. 
The fifth set of panels show that additionally excluding 
the [/-band lag from the fit shown in the fourth panel 
greatly reduces the xi although, as mentioned earlier, 
the fits are acceptable in all cases. This is discussed in 
Section 4.2. 

3.3. Spectral variability 

In this section we utilize the fact that emission from the 
AGN (central engine and surrounding regions) is vari¬ 
able while starlight from the underlying galaxy is not 
to separate these components and characterize the spec¬ 
tral shape of the AGN component. For most analyses, 


emission from the underlying galaxy is a complication to 
be removed from the UV/optical SFD before proceed¬ 
ing. One way to do this is image d ecomposition , as ha s 
been performed fo r NG C 5548 bv iBe ntz et al.l (j2009ll . 
iBentz et al.l (j2013| l and iMehdipour et al.l (j2015jl . Here 
we use an alternate approach to estimate the shape of the 
variable SFD (although not its normalization) directly 
from Swift data alone. 

We first filter the intensive monitoring data to include 
only orbits with observations in all six UVOT filters in 
order to obtain a uniform data set. For each band, we 
next measure the standard deviation of the flux (tr) and 
the mean error (e) and then calculate the error-corrected 
standard deviation, ac = \/cr^ — e^, although the error 
correction was always small, typically ^1%. This pro¬ 
vides a direct estimate of the intrinsic variability in each 
band but does not include the mean flux of the AGN. 
We then take the logarithm of both quantities and fit 
a function of the form logio(crc(A)) = alogiQ(A) -I- /3, 
which is equivalent to crc(A) oc A“, in order to measure 
the power-law slope a of the variable component. As 
shown in Figure 5, this fit yields a = —1.88 ± 0.20. In 
order to estimate the intrinsic shape of the variable AGN 
SFD, we perform the same exercise after first deredden¬ 
ing the data (assuming E{B — V)= 0.017; see Paper I). 
This yields a slope of a = —1.98±0.20. We note that this 
is consistent with predicted thin acc retion disk slopes of 
a = —2 to —2.33 (|Davis et al.ll2007ll . 

4. DISCUSSION 

4.1. Reverberation mapping of the accretion disk 

The standard model of a geometrically thin, optically 
thick AGN accretion disk predicts that the disk will be 
hotter in the inner radii and cooler in the outer radii, 
with dependencies on the black hole mass (and thus the 
Schwarzschild radius) and Fddington r atio. This is fo r 
instance quantified in Equation 3.20 of lPeters^ (1199711 . 

where T(r) is the temperature at radius r, M/M-Edd is 
the mass accretion rate divided by the Eddington rate, 
assuming a radiative efficiency of 77 = 0 . 1 , Mg is the 
black hole mass in units of IO^Mq, and i?sch is the 
Schwarzschild radius. 

Gombining Equation 1 with Wien’s law (Amax = 2.9 x 
lO^/T, where A is measured in Angstroms and T in 
Kelvin) and the idea that the lags are dominated by light 
travel times from the center (so r = r/c) yields the rela¬ 
tion T OC A^/^. Note that while this derivation was for a 
disk heated “internally” by viscous processes, the same 
X^/^ dependence will arise in a disk heated “externally” 
by the putative central corona . This is shown in Equa- 
ti ons 1 and 2 ofiGackett et al.l (1200711 and Equation 4.56 
of iNetzeil (|2013[ 1. However these derivations all assume 
a relatively flat disk. If the disk is stron gly warped or if 
the c orona is distributed across the disk (IDexter fc Agoll 
I 2 OIIII . then an externally-heated disk will be hotter at 
large radii, leading to a flatter lag-wavelength relation. 
The consistency of the observed r — A relation with the 
predicted relation broadly supports the standard accre- 
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tion disk temperature profile. 

Equation 1 can also be used to estimate source pa¬ 
rameters under the simple assumption of a face-on disk 
in which each annulus at temperature T{r) radiates all 
its luminosity at A^ax as given above. The observed 
relation r = — 1) has a single free pa¬ 

rameter B Ri 0.35 day, which indicates, for a face-on 
disk, that an annulus of radius ^0.35 It-day radiates at 
T = 2.9 X 10^/Ao = 2.2 x 10^ K for Aq = 1367 A. This 
distance of 0.35 It-day corresponds to r/i?sch = 90 for a 
3.2 X IO^Mq black hole. 

A more realistic picture would account for the fact 
that each annulus radiates as a blackbody of temper¬ 
ature T{r) instead of radiating at a single wavelength 
Amax- Accounting for that will yield a larger value for 
the radius at which the disk emission peaks at Amax be¬ 
cause more flux at a given wavelength is produced by 
hotter blackbody emission interior to radius R than by 
cooler blackbody radiation exterior to R. A future pa¬ 
per in this series (Starkey et ah, in prep.) utilizes direct 
modeling of each UV/optical light curve to produce a 
more rigorous analysis. 

Nonetheles s these large sizes m a y cau se problems for 
the standard iShaknra fc Snnva^ (jl973l l a-disk model. 
Assuming a value for M/M^dd = L/L-^dd = 0.03, de¬ 
rived using a disk luminosity L that is ^37% of the to¬ 
tal luminosity, yields an accretion disk radius of only 
r / Rsch = 40 using the same formula and assumptions as 
above. We note that there are models that produce ef¬ 
fectively larger disks which can potentially better explain 
the UV/optical variability pro perties of AGN, such a s the 
inhomogeneous disk model of lDexter fc Agoll (120111) . 

4.2. Contribution of BLR emission 

Another interesting result shown in Figure 4 is the 
longer {/-band lag, relative to the fit and to the lags 
of nearby bands. Excluding the {/-band data from the 
fit yielded a significant improvement in although 
the overall fit is acceptable in either case. The final 
fit on the bottom of Figure 4, which excludes {/-band, 
predicts a {/-band lag of r = 0.85 while the observed 
value is T = 1.35 ± 0.24, a difference of 2 ct. In retro¬ 
spect, one can see in both th e previous NGC 5548 cam¬ 
paign ( McHajd^_et_jMj201^ and the NGC 2617 cam¬ 
paign ( Shapoee et al. 1 120141 1 that the {/-band lags were 
larger than the R-band lags, although those campaigns 
measured lags with much larger errors, so the deviation 
was not significant. The vastly superior short timescale 
sampling provided by the current campaign (see Fig¬ 
ure 1) allows for the measurement of this apparent effect 
with higher significance. 

There is a simple explana t ion fo r this excess lag, dis¬ 
cussed by iKorista fc GoadI (|20nih : Balmer continuum 
emission (both thermal diffuse and reflected incident con¬ 
tinuum) and other pseudo-continuum emission from BLR 
clouds (e.g., UV Fe ii) contributes significantly to the 
observed {/-band flux, and since the BLR is much larger 
than the optically-bright accretion disk, it will increase 
the ob served lag. This effect was seen by iMaoz et al.l 
(|1993ll in NGC 5548 RM data from the 1989 campaign. 
The sensitivity of the strength of the diffuse continuum 
component to the presence of high gas densities and high 
ionizing photon fluxes make it an important diagnostic of 


the physical conditions within the BLR. A possible alter¬ 
native is that the Balmer continuum (and other pseudo- 
continua) is produced in an “intermediate” region smaller 
in size than the classical BLR but larger than the accre¬ 
tion disk. This could more naturally explain the rel¬ 
atively small increment in the {/-band lag, although it 
would also mean adding a previously-unknown emission 
component to the many already required to explain AGN 
spectral energy distributions. 

We note that continuum light curves measured at 
longer wavelengths and/or with narrower bands will be 
much less sensitive to this effect. This cannot be done 
with Swift, but future papers in this series will ana¬ 
lyze an expanded set of HST and optical photometric 
bands (Fausnaugh et al. in prep.) and ground-based 
spectroscopy of NGC 5548 (Pei et al. in prep.), providing 
a more sensitive test of the degree to which these CCFs 
are contaminated by emission from hot gas surrounding 
the central engine. 

4.3. The aecretion disk size - black hole mass relation 

As discussed in Section 1, quasar gravitational mi- 
crolensing studies have been used to estimate accre¬ 
tion disk sizes, finding a tenden cy for disk sizes to 
increase with black hole mass (jMorgan et al. I 120101 
IMosauera et al.l 1201311 . Disk RM measurements of 
Seyfert galaxies can be used to extend such relations 
to lower masses, luminosities and (probably) Edding¬ 
ton ratios, generally with smaller uncertainties because 
of the greatly reduced physical complexity of the mea¬ 
surement. Figure 6 shows a summary of microlensing 
sizes estimates (the half-light radius i?i /2 at rest frame 

2500A) from IMosauera et al.l (|2013ll as open triangles. 
Microlensing studies frequently focus on R1/2 because es¬ 
timates of its value are relatively insens itive to changes in 
the u nderlying (disk) emission profile (IMortonson et alH 
1200511 . A fit to these data as a power-law, logiQ(i?i/ 2 ) = 
A -I- BlogiQ(M/Mo), with Mg = 3 x IO^Mq to mini¬ 
mize covariances between the parameter estimates and 
assuming 0.3 dex uncertainties in the black hole mass 
estimates, yields A = log(i?o/cm) = 15.81 ± 0.16 and 
B = 1.29 ± 0.33 with = 8.73 for 9 dof, a statistically 
acceptable fit. The slope of the fit is driven by the higher 
mass systems, leading it to lie below the measurements 
in the mass range of NGC 5548. 

The comparable disk RM size from the present study 
is the distance corresponding to the lag for the UVWl 
filter centered at 2600A, which we show as the filled 
square labeled “N5548”. This combines the parameter B 
(= 0.35 ± 0.04 days for the simple face-on model in Fig¬ 
ure 4) with the HST UV to UVWl lag (0.61±0.20 days), 
to give an estimate that i?i /2 = 0.96 ± 0.21 light-days. 
This assumes that the size corresponding to the observed 
lag corresponds to i?i/ 2 , which may not be correct, but 
is a reasonable assumption pending a theoretical model 
for how the disk RM lag should be interpreted in detail. 
For example, R1/2 is 2.44 times larger than the radius at 
which the photon wavelength matches the disk temper¬ 
ature discussed in Section 4.1. A similar procedure was 
used to add the UVW l size estimate for NGC 2617 from 
iShappee et al.l (|2014f) , where the error bar is designed to 
span their systematic uncertainties. 

While there are residual systematic uncertainties in 
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this comparison, such as the meaning of the disk RM 
lag as a physical size and potential differences in the Ed¬ 
dington ratios of the nearby lower luminosity Seyfert Is 
and the distant high luminosity quasars, the results from 
the two very different methods are broadly consistent. If 
we simply fit the combined data, we find A = 15.96±0.12 
and B — 0.98 ± 0.23 with = 13.51 for 11 dof, which 
is shown by the dashed line in Figure 6. The slope is 
flatter and better defined, the fit is consistent with the 
combined data, and the parameters are consistent with 
the results using only the microlensing results. Inter¬ 
estingly, the slope B is almost exactly unity (and thus 
the fit line is nearly parallel with the last stable orbit, 
shown as a solid line in Figure 6), indicating that the 
disk size in units of Schwarzschild radii is nearly con¬ 
stant, i?i/2(2500A) ^ lOORsch, over a very wide range 
of AGN masses. That these two radically different meth¬ 
ods agree this well seems remarkable given the different 
underlying physics and the possible range of systematic 
effects. 

4.4. Relation of X-ray to UV/optical continua 

Figure 3 shows that compared to the strong (r^ax = 
0.57— 0.9) correlations within the UV/optical, the cor¬ 
relation between the HST 1367 A and the X-ray bands 
is much weaker (rmax = 0.35 — 0.44). This is surpris¬ 
ing because it is well established that the optical and 
X-ray light curves of NGC 5548 are very w ell correlated 
(rmax = 0.95) on longer timescales of years (|Uttlev et al.l 
1200311 . That is, the strong long timescale optical/X-ray 
correlation does not translate to strong short timescale 
UV/X-ray correlations in NGG 5548. Visual examina¬ 
tion of the NG C 5548 UV and X-ray li ght curves both in 
this paper and iMcHardv et al.l (|2014[1 shows that there 
are some periods in which the UV appears to lead the 
X-rays, some in which the UV appears to lag the X-rays, 
and some in which there is no simple discernible relation¬ 
ship. 

Periods of uncorrelated X-ray/optical v ariations are 
also s een in other Seyferts (e.g. NGC 3516. lMaoz et al.l 
120021 and Mkn 79. lBreedt et al.ll2009| l. The phenomenon 
may be linked to internal heating fluctuations in the disk 
which are not ‘seen’ by the X-ray emitting region, per¬ 
haps linked to mass accretion fluctuations which do not 
propagate to the central X-ray emiting region due to 
viscous damping (e.g., the explanation of a similar phe¬ 
nomenon seen on equivalent, mass-scal ed time-scales in 
a ste llar mass black hole X-ray binary iCassatella et al.l 
l2012f) . This may indicate that a significant fraction of the 
UV/optical emission is not due to reprocessing of X-ray 
photons, but rather is generated internally. In this case, 
the observed time lags would not be dominated by light 
travel effects, but instead would depend on the physics 
of the internal disk variations. 

Alternatively, the lack of correlation between the 
UV/optical and X-ray could be a signature of absorp¬ 
tion due to intervening material, but then one would 
expect the UV to be better correlated with the hard X- 
rays than the soft X-rays, because the latter would be 
much more strongly affected (and the light curves more 
decorrelated) by ’’warm” absorption. We do know that 
NGC 5548 shows strong var iable absorption in the X- 
rays (|MehdiDOur et al.ll2015ll . and there are indications 


that the absorption was changing (decreasing) during the 
Swift campaign. Nonetheless the fact that the hard X- 
rays show a smaller correlation coefficient than the soft 
X-rays suggests that this may not be a complete explana¬ 
tion. A somewhat different scenario, based on a correla¬ 
tion between_soft2Gray_excess and far-UV also observed 
by (iMehdipour et al.ll2015ll . is that both are associated 
with Comptonization. Finally, it could simply mean that 
the observed 0.3-10 keV X-ray band is a poor proxy for 
emission from the putative hot corona, which should emit 
the bulk of its luminosity at harder energies. At this 
point it is not possible to say with certainty which, if 
any, of these explanations is correct. 

5. CONCLUSIONS 

This paper presents the results of the most intensive X- 
ray/UV/optical AGN monitoring ever, spanning a dura¬ 
tion of months. We find that the UV/optical light curves 
are all well correlated with lags of ~l-2 days increasing 
to longer wavelengths. These lags are well-fitted by the 
relation r oc in agreement with standard steady- 

state accretion disk predictions under the assumption 
that time lags are dominated by light travel times. The 
fits yield a disk size of ^ 0.35 ± 0.05 It-day at 1367 A, 
larger than expected from standard a-disk models or ex¬ 
trapolation from higher-mass microlensing studies. In¬ 
terestingly the [/-band lag is anomalously large, suggest¬ 
ing that the U band is affected by Balmer continuum 
emission from the BLR. The X-ray/UV correlations are 
weaker and less consistent, however, so these data do not 
confirm all predictions of the reprocessing picture. 

We are planning a series of future papers to explore 
these results in greater detail. Fausnaugh et al. (in prep.) 
will present ground-based optical and further HST con¬ 
tinuum data, allowing a check on the wavelength depen¬ 
dence of the observed interband lags. Pei et al. (in prep.) 
will use ground-based spectroscopy to measure the con¬ 
tinuum in narrow spectral windows much less affected 
by BLR emission, further refining this analysis. Starkey 
et al. (in prep.) will apply Markov Chain Monte Carlo 
methods to directly model disk emission from these con¬ 
tinuum data, allowing a much more direct probe of the 
physical conditions. 

This RM disk size estimate of a relatively low-mass 
Seyfert 1 galaxy forms a nice complement to more nu¬ 
merous but more uncertain accretion disk size estimates 
derived from generally higher-mass quasar microlensing 
studies. The combination of the two datasets allows im¬ 
proved determination of the accretion disk size - black 
hole mass relation, which interestingly suggests 2500 A 
accretion disk sizes of i?i /2 ~ lOORsch- Further Seyfert 1 
accretion disk RM experiments will allow this relation to 
be tested and refined. 

Most important for the long term is that this exper¬ 
iment demonstrates how dense and broad coverage in 
both wavelength and time can be used to probe a nearby 
AGN accretion disk with unprecedented detail. Swift was 
named “the premier facility for multi-wavelength time 
domain astronomy” by the latest NASA Senior Review 
PanelI3 There is certainly no other observatory that can 
single-handedly monitor AGN with such dense and broad 

http://science.nasa.gOv/media/medialibrary/2014/05/15/ 
Final_Report_Astro2014_SeniorReview_Panel.pdf 
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temporal and frequency coverage in the UV/optical and 
X-rays. We expect that this experiment will become a 
template for future Swift campaigns that characterize the 
accretion disks of a sample of AGN covering a range of 
black hole masses, Eddington ratios and other source pa¬ 
rameters. 
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APPENDIX 

As discussed in Section 2.2 we discovered “dropouts” in the UVOT light curve in the course of the data reduction: 
isolated points with fluxes many sigma below those of their nearest neighbors. Figure A1 shows the UVOT light curves 
of NGC 5548 after initial flux measurements and removal of data points affected by tracking problems. The dropouts 
are the points with red error bars, most frequently seen in the UV bands. In order to quantify this effect, we first 
parametrized the deviation for every point in the light curve as df = (#v — 0.5 x (Fat-i -I- where F/v 

and (Tat are the measured flux and Icr error bar for visit N. To minimize the effects of intrinsic variability, we only 
tested data with observing gaps tN+i — tx-i < 2.5 days is the THJD of the iVth visit); all visits with longer gaps 
between their nearest neighbors were ignored. 

We then flagged all points with negative excursions greater than the largest positive excursion seen in that filter as 
dropouts. The largest positive excursions used to define the threshold of what is a dropout are themselves sensitive to 
dropout measurements, in that the largest df values tend to be found when the TV — I or V -|- I flux measurement is 
a dropout. We therefore redefine these thresholds iteratively, removing the dropout points and then re-evaluating the 
largest positive excursions. This is repeated until the largest positive excursion remained the same, that is, it was not 
associated with a dropout. This procedure limits the number of false positives to of order one per light curve. Note 
the strong dependence on UVOT band, with 33, 13, 21, II, 2 and 1 dropouts for the UVW2, UVM2, UVWl, U, B 
and U bands, respectively. 

We then mapped the source location of every exposure for the three UV filters {UVW2, UVM2, UVWl) to the 
UVOT detector coordinates (Figure A2). Blue dots show the points used in the analysis, red Xs the dropouts as 
defined above, and open black circles the points ignored because they lacked nearby (in time) neighbors. Figure A3 is 
a blow-up of the region with the dropouts. Note that the dropouts cluster together such that most can be enclosed 
by a small number of boxes. We define eight rectangles in the detector plane, one for each cluster of at least three 
dropout points. The coordinates of these boxes are given in Table Al. We then went back to all six filters (including 
the optical UBV bands) and flagged every visit/filter that falls in any of these boxes. 

These additional suspect points are shown as red Xs in Figure Al. We ran Kolmogorov-Smirnov tests on each band 
to test if the deviation for points in the suspect regions (the red Xs in Figure Al) derives from the same population 
as the unaffected data (black dots in Figure Al). We find that for five of the six UVOT filters (all except V) the two 
samples are not consistent at the 10“"^ to 10“^® level, with the strongest differences at the highest frequencies. We 
then eliminated all data points in any filter that fell in these eight boxes. These points are shown as red Xs in Figure 
Al. The remaining points (the black dots in Figure Al) then formed the final light curve shown in Figure 2 of the 
main section of this paper. 
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TABLE Al 

UVOT DETECTOR BAD BOXES 


xl 

x2 

yi 

y2 

357 

371 

634 

651 

415 

440 

582 

632 

441 

469 

634 

655 

448 

457 

440 

442 

479 

484 

539 

542 

516 

537 

596 

619 

536 

538 

564 

568 

545 

583 

573 

607 


Note. — UVOT detector coordinates 
(xl,x2,yl,y2) for the eight “bad” boxes 
shown in Figure A3. These are image co¬ 
ordinates for a full-frame raw UVOT im¬ 
age with 2x binning (the default UVOT 
image mode). 
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Histogram of UVW2/U pairs 
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Fig. 1.— Histograms showin g the number of UVW2 /U pairs for the current intensive NGC 5548 monitoring cam paign (top), the earlier 
NGC 5548 campaign ('middle. IMcHardv et aLll2014l) and the NGC 2617 campaign (bottom, IShappee et al. 1120141) . Data are binned by 
orbit, and all pairs with separations of less than half an orbit are excluded. UVW2 was used because it was the most frequently observed 
UVOT band, while U, a typical less well sampled band, was used because that band is particularly interesting (see Section 4.2). The range 
±2 days is shown because this is the key cadence range that has not previously been well sampled. Note that the current campaign samples 
these short cadences ^^5-8 times more frequently than the previous NGC 5548 and NGC 2617 campaigns. 
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Fig. 2.— Light curves for the intensive monitoring period (HJD 2,456,706-2,456,831), going from shortest wavelength (top) to longest 
(bottom). The band name and central wavelength are given on the left of each panel. Top two panels show the Swift hard and soft X-ray 
(HX and SX respectively) light curves, in units of c/s. Third panel shows the HST light curve, in units of 10~^‘^ergs“^cm“^A“^. Error 
bars for this light curve are typically '^1.5%, just barely visible in the plot. The bottom six panels show the Swift light curves, again in 
units of 10“^'^ erg cm“^s“^A“^ . Dashed gray lines show times THJD 747.179, 785.752 and 818.993, three local maxima of the HST light 


curve. 
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Fig. 3.— (3a) Interpolated cross-correlation functions for the intensive monitoring period light curves (Figure 2), with all correlations 
measured relative to the HST light curve, after removing long term trending (see Section 3). The band name and central wavelength 
are given on the left of each panel. Note that the interband lag goes from negative to increasingly positive as the band’s wavelength 
increases. Note also that the UV/optical correlations are all strong {rmax = 0.57 — 0.90) but the X-ray/UV correlations are much weaker, 
(Xmax < 0.45). (3b) Cross-correlation centroid histograms derived from the CCFs as discussed in the text. The band name and central 
wavelength are given on the left of each panel. All distributions except HX appear consistent with a Gaussian. 



































14 


Edelson et al. 






Fig. 4.— Lag-wavelength fits based on the data in Tables 1 and 4. The top row shows the most general fit, r = A-\- — 1), with 

the power-law index C allowed to float. The next row fixes the index at the theoretically expected value C = 4/3. All data are included in 
the first two sets of fits. In the third row the intercept is fixed at A = 0 and the HST ACF data are excluded. The fourth row shows these 
fits with the X-ray data HX and SX also excluded, while the fifth row additionally excludes the f/-band data. See text for further details. 
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Fig. 5.— The error-subtracted variable flux (cr(; 7 (A)) as a function of wavelength (A). The original data are shown as red Xs and the 
dereddened data as blue crosses. A fit to the function logio(f 7 c(A)) = Q:logj^Q(A) -|- /3 yielded a = —1.88 ih 0.20 for the original data (solid 
red line) and a = —1.98 ib 0.20 for the dereddened data (dashed blue line). 



Fig. 6. — Accretion disk size estimates from quasar microlen sing studies (open tr iangles, IMosauera et al.11201311 . the current study of 
NGC 5548 (black filled square) and NGC 2617 (open square, IShapoee et al.f|2Q14l') as a function of black hole mass. The dotted line 
shows the fit to just the microlensing data and the dashed line the fit to all data points including the Seyfert 1 RM measurements. To 
give a sense of other scales associated with accretion disks, the lower solid line shows the last stable orbit of a non-rotating black hole at 
R = 3-Rsch = QGMjc^ and the upper solid lines shows R 1/2 at 2500A for a simple thin disk with Eddington ratio = 0.03 and 

efficiency 77 = 0.1 to provide a sense of scale. The latter curve can be shifted as {L/L-Edd)^^^ tor different choices of these factors. 



















16 


Edelson et al. 



400 500 600 700 800 

Truncated Heliocentric Julian Day Number (THJD = HJD - 2,456,000) 


Fig. Al.— Initial UVOT light curves of NGC 5548 for the period THJD 390-832. The error bar colors indicate the results of the dropout 
test: black errors indicated that the point passed the dropout test (small deviation), red error bars indicated that it failed the test (large 
deviation) and cyan error bars indicate that it was not tested (as it lacked sufficient nearby neighbors). The symbols show if the point fell 
inside/outside the UVOT boxes (shown in Figure A3): black dots fell outside the boxes and thus were used in the final light curve (shown 
in Figure 1) and red Xs fell inside the boxes and were excluded from the final light curves. 
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Fig. A2.— UVOT detector coordinates of the UVW2, UVM2 and UVWl data in Figure Al. Data that were tested for dropouts are 
shown as blue dots and those that were not tested (due to their not having sufficiently nearby neighbors) are shown as open black circles. 
The points that failed the dropout test are marked with red Xs. The black rectangle denotes the region shown in Figure A3. 
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Fig. A3.— The detector region containing all the observed UVOT dropouts. The black circles, blue dots and red Xs have the same 
meaning as in Figure A2. Green rectangles are drawn around clusters of three or more dropouts. Eight such rectangles enclose most red 
Xs and are used to define suspect regions on the UVOT detector. All data points within these regions were eliminated to form the final 
light curve (Figure 2), regardless of whether or not they failed the dropout test. 

































































































